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We investigate the effect of the electric-charge neutrality in /3-equilibrium on the chiral phase 
transition by solving the chiral and diquark condensates in the two-flavor Nambu-Jona-Lasinio 
model. We demonstrate that the electric-charge neutrality plays a similar role as the repulsive 
vector interaction; they both weaken the first-order chiral phase transition in the high-density and 
low-temperature region. The first-order chiral phase transition is not affected, however, at finite 
temperatures where the diquark condensate melts. In this way the chiral phase transition could 
be second-order at intermediate temperatures if the diquark effects overwhelm the chiral dynamics, 
while the first-order transition may survive at lower and higher temperatures. The number of the 
)sL} critical points appearing on the phase diagram can vary from zero to three, which depends on the 

, relative strength of the chiral and diquark couplings. We systematically study the possibility of 

the phase structure with multiple QCD critical points and evaluate the Meissner screening mass to 
confirm that our conclusion is not overturned by chromomagnetic instability. 
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, Cti . PACS numbers: 12.38.Aw, ll.10.Wx, 11.30.Rd, 12.38.Gc 

m 

It is generally believed that quantum chromodynamics (QCD) exhibits a rich phase structure in extreme environment 
^"'i such as high temperature and high baryon density. In the last decade the color-superconducting (CSC) phase has 
attracted lots of theoretical interests and triggered extensive studies of dense and cold quark matter p| 0, [H, 01. At 
asymptotically high density that justifies the perturbative QCD calculations the color-flavor locked (CFL) phase Q has 
been established as the ground state of quark matter. There, the pairing energy, the critical temperature, the screening 
properties, the collective excitation energy, and so on are well understood from the first-principle calculations. 

In a practical sense, however, the accessible baryon density in nature would be, at most, ten times the normal 
| nuclear density even in the interior of the compact stellar objects. In terms of the quark chemical potential the 
corresponding value should be less than 500 MeV in reality and thus other MeV energy scales such as the current and 
dynamical quark masses, the electric chemical potential, etc take part into the dynamics. We are yet far from thorough 
understanding on the phase structure in this density region which is commonly referred to as the intermediate density. 
Going down from the CFL phase, we have to rely on chiral effective theories to unveil the possible CSC phases and 
OO , to draw the boundary lines on the phase diagram which separate different CSC states [1,0). 

Recently, an interesting proposal has been made for the plausible phase structure in the intermediate density 
region Q. That is, there may exist a new QCD critical point (i.e. the terminal of the first-order phase boundary) 
induced by the UA(l)-breaking vertex. This speculation is based on a general Ginzburg-Landau theory in terms 
of the order parameter fields constrained by QCD symmetries. There, the three-flavor anomaly term generates the 
coupling between (^ip) (chiral condensate), (ipif)) (diquark condensate), and ($r/j) (antidiquark condensate). It has 
been then argued that the resultant crossover of chiral restoration at small temperature embodies the hadron-quark 
continuity hypothesis Q • Whether the new critical point appears with reasonable parameter set or not needs further 
investigation. Maybe, to settle the situation without ambiguity, we should wait for future developments in lattice 
QCD simulations at finite density. 

In a different context, before discussed in Ref. Q, the appearance of another critical point at low temperature was 
pointed out in Ref. [1] within the two-flavor Nambu-Jona-Lasinio (NJL) model. The crucial ingredient in Ref. Q is 
the four-fcrmion interaction in the vector channel as well as in the scalar one. In this case, the vector interaction 
diminishes the first-order chiral phase transition [{J [l(| [H| , which allows for enhanced competition between the chiral 
and diquark condensates in the widened coexisting phase. For some choices of the NJL model parameter two critical 
points show up along the phase boundary that signifies crossover, first-order transition, and crossover again with 
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increasing temperature. We note that a similar phase diagram with two critical points is suggested in the lattice 
calculation of two-color QCD [13], though it seems to have not been established yet. As mentioned in Ref. Q, the 
vector-channel interaction could give a possible explanation for this two-color phase structure with two critical points. 

In this paper we will reveal another mechanism leading to multiple QCD critical points; the electric-charge neutrality 
realized in (3 equilibrium can make chiral restoration being smooth crossover at low temperature, so that an analogous 
situation to Refs. 0, [1[ takes place. It is worth remarking here that imposing the electric neutrality has a similar 
effect on the thermodynamic potential to introducing the vector-channel interaction. One can understand this in 
terms of the electric chemical potential, p e , which we define such that p e > not for positron but electron. In the 
model treatment without gauge fields a nonzero p e mimics the role of Aq to neutralize the system. The point is 
that p e enters the dynamics like the vector-channel interaction; a vector-channel interaction — Gy('07 At V') 2 induces an 
effective (renormalized) chemical potential [pj which takes a form, 

/!R = ih - 2G V • p q ■ (1) 

Here Gy > represents the repulsive vector coupling constant and p q is the quark number density, (ipj ip). For 
two-flavor QCD p, e leads to a mismatch between the M-quark chemical potential p u = p — |/i e and the d-quark one 
P-d = p + \p e - We shall use a notation, /2, to denote the average chemical potential for it and d quarks resulting in 

P = P - \pe ■ (2) 
6 

Comparing the forms of /zr and p above, one may well anticipate that p e ^ can be taken in effect as a repulsive 
vector coupling for the bulk properties, in addition to keeping the electric neutrality. In physics terms, the electric 
chemical potential realizing a finite electric-charge density in part plays a role as a baryon chemical potential on its 
own. This situation is drastically different from the massless three-flavor case where the electric charge generator 
happens to be traceless so that no coupling between the electric-charge and baryon fluctuations arises (i.e. no term 
proportional to p e arises in p in this case). 

In the best of our knowledge the direct coupling between the quark density and the electric chemical potential 
has drawn only little attention so far, though the vector-channel interaction and the electric neutrality have been 
investigated in separate contexts. In view of the results in Ref. Q, which is produced solely by the vector interaction, 
it is natural to expect that p e may also have a significant impact on both the chiral and CSC phase transitions. 

The purpose of this paper is to investigate this issue seriously and depict an intuitive picture which opens a 
possibility to drive more QCD critical points than only one. We shall demonstrate our idea in the framework of the 
two-flavor NJL model for concreteness. Indeed, the competition between the chiral and diquark dynamics results in 
zero, one, two, and three critical points depending on the relative strength of the chiral and diquark couplings. 

The paper is organized as follows. In Sec. [HI the model is introduced and the formalism is presented. The numerical 
results and discussions are given in Sec. IIIII The final section is devoted to the summary and concluding remarks. 

II. MODEL AND FORMALISM 

We will explain the choice of the effective model, the model parameters, and the resulting thermodynamic potential 
in order. Since we adopt a standard description by the NJL model, the experts could skip to our results in Sec. IIIII 

A. Model 

Varieties of NJL-type models have been extensively used to investigate the CSC phase transition at moderate and 
large density as well as at zero density [H, [3]- For the two-flavor case, the commonly used Lagrangian of the NJL 
model reads 

£ N JL = i> (*7 M ^ ~ ™0) lp + £qq+ C qq , (3) 

where the chiral interaction part is 

£ qq = G[(^) 2 + (^ 5 ^) 2 ] , (4) 

and the diquark part which is relevant to the mean-field condensate (i.e. the spin, flavor, and color are all antisym- 
metric) is 

C qq = H [(^C 15 t 2 X a ^ t )(^ t C 75 t 2 X a ^)] . (5) 
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Here, C — 17072 stands for the Dirac charge conjugation matrix, and G and H are the coupling constants for the 
mesonic and diquark channels. The current quark mass matrix is given by m = diag(m u , m c i) in two flavors and we 
shall work in the isospin symmetric limit with m u — rrid — m. We note that A^'s are the antisymmetric Gell-Mann 
matrices (i.e. A runs over 2, 5, 7 only) for the color SU(3) group and •Ps are the Pauli matrices in flavor space. 

For simplicity, the coupling constant in the vector channel is set to zero in this paper. The effect of the electric- 
charge neutrality on the chiral phase transition with including the Polyakov loop dynamics and the nonzero vector 
interaction will be reported in our future work fl5T ] . 

It deserves notin g he re that the scalar four-fermion interaction Cq q in general consists of two types of different 
interactions [l3l [l4l.ll6l]. that is, 

d = Gi [(^) 2 + (V^V>) 2 + fe^) 2 + (^75T<A) 2 ] , (6) 
C 2 = G 2 [(^V) 2 ~ WN>) 2 - (V^75V>) 2 + i^l^f] ■ (7) 

Both interaction terms have the symmetry of SU(2)l X SU(2)r, x U(l), while the axial symmetry, Ua(1), remains 
only in C\. The Ua(1) breaking part C 2 belongs to the instanton-induced (two-flavor 't Hooft) interaction. In the 
mean-field level the constituent quark masses in the presence of L\ and C 2 are 

M u = m-4G 1 (uu) - 4G 2 (dd) , M d = m - 4G?i (dd) - 4G 2 (uu) . (8) 

Therefore, in general, M u ^= Aid if there exists a chemical potential mismatch between u quarks and d quarks by /i e 
and thus (uu) 7^ (dd) . If we introduce a parameter a to relate G\ and G2 to G in a way that 

Gi = (l-a)G, G 2 = aG, (9) 

we notice that the standard Lagrangian ([4| corresponds to the case of a — 0.5. In such a case, as is clear from Eq. (|5J), 
the constituent mass of u quarks is always identical to that of d quarks regardless of a difference in fi u and /id- Once 
we get ready to proceed to the numerical calculations, in Sec. Ill CI we will check the dependence on a in a simple 
case without diquark condensation. In any case, because this paper aims to illustrate a general feature in the phase 
structure, we shall stick to the simplest choice a = 0.5 which makes no difference in the qualitative picture. 

Then, there are four model parameters left; the current quark mass m of u and d quarks, the coupling constants G 
and H, and the three-momentum cutoff A. In this work, we take the same parameters as in Ref. |17j | which are fixed 
so as to reproduce the three physical quantities in vacuum; the pion mass ~ 140 MeV, the pion decay constant 
/tt w 94 MeV, and the chiral condensate (uu) = (dd) w -(251 MeV) 3 with 

m = 5.5MeV, G = 5.04GeV~ 2 , A = 0.651GeV. (10) 

The corresponding constituent quark mass in vacuum is 325.5 MeV for this set of the model parameter. The standard 
value of the ratio H/G is 3/4 = 0.75, which is deduced by the Fierz transformation from the local current-current 
interaction. In this paper, we rather treat this ratio as a free parameter and shall perform a systematic survey. 



B. Thermodynamic Potential with Neutrality Condition 

In general, the quark chemical potential matrix jj, takes the form [l8T | 

p, = n - [i e Q + jti 3 T 3 + ^ 8 T 8 , (II) 

where [i is the quark chemical potential (i.e. one third of the baryon chemical potential), \x e is the chemical potential 
associated with the (negative) electric charge, and ^3 and /ig represent the color chemical potentials corresponding 
to the Cartan subalgebra in color SU(3) space. The explicit form of the electric-charge matrix is Q — diag(|,— |) 
in flavor space, and the color charge matrices are T 3 = diag(i, —^,0) and Tg = diag(-|, 5, — |) in color space. The 
chemical potentials for different quarks are listed below; 

211 211 

Mru = M — jjMe + 2^3 + 3^8 j Mgu = M ~ 3/^ _ 2 ^3 + 3M8 j 

Urd = /-!+ \u e + |/i 3 + , Vgd = M + jMe - 5M3 + 5/^8 , (12) 

2 2 12 

Mim — I 1 — 3^6 — 3M8 j Mbd = M + ijMe — 3M8 ■ 

The four-quark interactions develop a dynamical quark mass with nonzero chiral condensate as 

M = m-a = m- 2G(^) , (13) 
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while the diquark condensate A and antidiquark condensate A* could appear at high enough baryon density. Here, 
we follow the common treatment for two-flavor CSC that the blue quarks do not take part in the Cooper pairing. 

Using the standard bosonization technique, the mean-field thermodynamic potential in the NJL model with the 
diquark degrees of freedom as well as the electron contribution takes the following form: 

fi = — + — L_ (n 4 + 2n 2 T 2 u 2 + — T 4> ) -TV [ ^ TV In ^A^gj ( U ) 

4G + AH 12n 2 + Zn ^ + 15 J {2nf T ' [ ' 

where the sum runs over the Matsubara frequency lo„ — (2n + l)nT and Tr is taken over color, flavor, and Dirac 
indices. The inverse quark propagator matrix including both the chiral and diquark condensates in the Nambu-Gor'kov 
formalism is then given by 

= ( i^;' A2 t<v-' 2 1 ■ (15) 

[Go]' 1 =7o(«^n±A) -J-p-fh. (16) 
Taking the Matsubara sum, we can express the thermodynamic potential as usual as 
il(fj, e ,fi 3 ,fi & ,a, A;/i,T) 

* ' ^ + 2 ^ + ^_y;/' 7 ^ {Ei + 2Tln(1 + e -^ )) , <"> 



with 



4G AH 12tt 2 V^ e ^ e 15 



^)-g/(^< a + 2rl " <1 



with the dispersion relations for six quasiparticles [that is, 2 flavors x 3 colors; the spin degeneracy is already taken 
into account in Eq. (|17[) ] and 6 quasianti-particles. The unpaired blue quarks have the following four energy dispersion 
relations, 

Ebu = E — (ibu , Efyu = E + fj.hu Ebd = E — [ibd ,Ebd~E + fibd (18) 
with E = \J p 2 + M 2 . In the rd-gu quark sector with pairing we can find the four dispersion relations, 



E rd-gu =E A ± §(/Vd ~ Hgu) = E A ± |(/i e + /i 3 ) : 
E rd-gu =E A ± \{^rd ~ Hgu) = E A ± |(/i e + /i 3 ) . 

and the ru-gd sector has another four as 

E ru-gd = E A ± \ {p ru - fXgd) = E A =p |(/i e ~ Ms) • 



(19) 



(20) 



where _Ea = y/(E — n) 2 + A 2 and i?A = y/(E + /J.) 2 + A 2 . The average chemical potential is defined by 

^ = ^^ = ^^ = /i -y + y (21) 

For the two-flavor CSC case, the color charge corresponding to the matrix T3 is always zero since the color SU(2) 
symmetry is left unbroken for red and green quarks. That means /13 = 0. In contrast to that in the NJL model, 
nontrivial coupling to the Polyakov loop might induce a nonzero /13 in the case of the Polyakov loop augmented NJL 
(PNJL) model [l9|,|20|], which is beyond the current scope. 

Since we know that /is is much smaller than fi e [l8l. |2~H to neutralize two-flavor CSC matter, the positive /i e /6 is 
overwhelming in Eq. (|2ip so that we can neglect fj,g in the numerical calculation. The average chemical potential then 
amounts to Eq. |(3J|. 

Minimizing the thermodynamic potential (|17|) . we can solve the mean fields a and A together with the chemical 
potential /i e from 

m = on = an 

da OA dfJ, e [ ' 
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FIG. 1: The constituent quark mass difference between u and d quarks for a = 0, 0.1, 0.2, 0.7, 0.9, where a is a parameter 
to indicate the flavor- mixing interaction defined in Eq. @. All results are obtained under the electric charge neutrality and 
without diquark condensation. 



150 



120 



> 

CD 




260 



280 300 320 340 360 
jlx [ MeV ] 



FIG. 2: The phase diagram for the chiral phase transition under the electric-charge neutrality. The solid line (dashed line) 
represents the first-order transition (smooth crossover) and the filled circle dot locates the QCD critical point. 



C. Dependence on a and the Constituent Quark Mass Difference 



For the cases with a 7- 0.5, as we have mentioned, M u should be different from Md in the presence of nonzero /z e , 
which is apparent in Eq. ©. Figure Q] shows the dependence on a in the behavior of M u , Md, and fi e as a function 
of \i at a fixed temperature T — 5 MeV. We do not take account of the diquark condensate for the moment. From 
the figure we note the following two points: First, ^ e becomes larger with increasing a, and the first-order phase 
transition tends to occur at a higher chemical potential. This means that the value of /i e , which plays a role similar 
to the vector interaction in the chiral phase transition, is sensitive to the magnitude of the flavor-mixing interaction. 
Second, in contrast to their absolute values of masses, we see that the mass differences between u and d quarks are 
sizable for a — and a — 0.1, while the mass difference becomes minor for a > 0.2 . 

It should be noted that the strength of the flavor- mixing interaction which originates from the Ua(1) anomaly may 
be small in comparison with the U(2) x U(2) symmetric part, and it would be intriguing to explore its effect on M M and 
Md- We stress, however, that we should perform such studies for the three-flavor case where the quantitative effect 
of the Ua(1) anomaly is more clearly seen in the rj-rj' system (for instance 22]). In the present two-flavor analysis, 
therefore, we shall only consider a = 0.5 for elucidating the effect of /i e on the chiral and CSC phase transitions. 
Actually, Fig.[T]shows that we can reasonably ignore the difference between M u and Md induced by fi e unless a < 0.2. 

It is known that as the coupling constant Gy of the vector interaction increases, the whole critical or crossover 
line of the chiral phase transition shifts toward larger chemical potential and the critical point moves toward smaller 
temperature and larger chemical potential, which disappears eventually at a large value of Gy [j| [l(| [HI- Let us 
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FIG. 3: The schematic change of the phase structure with increasing H/G from the left to the right. The solid and dashed 
lines represent the first-order transition and crossover, respectively. The number of the critical points depends on H/G. 

show that the fi e affects the phase diagram in the similar way as Gy, which is anticipated from Eqs. (fTJ) and $2$. To 
demonstrate it clearly, we try to enhance the effect of \i e artificially by varying the it-quark electric charge Q u by 
hand as Q u — |e|, 2|e|/3, |e|/3, while keeping the d quarks unchanged: The case of Q u = |e|/3 corresponds to the 
situation with no net effect on the average chemical potential. The real world is Q u — 2|e|/3. We can induce a further 
large net shift in the average chemical potential by taking Q u — \e\. Figure [2] shows that the critical point shifts 
from (T» = (36 MeV, 332 MeV) for Q u = |e|/3 to (T,/z) = (24 MeV, 345 MeV) for Q u = 2|e|/3, and eventually 
disappears from the phase diagram when we choose Q u — |e|. These results are quite reminiscent of the effects of 
changing Gy discussed in the literature [1, 0, [l(| El • 

III. NUMERICAL RESULTS AND DISCUSSIONS 

In this section, we shall discuss the effect of electric-charge neutrality in /3-equilibrium on the phase structure with 
the CSC phase as well as with the chiral transition taken into consideration. We shall show that the /j e induced by 
the neutrality constraint gives rise to a phase structure with multiple critical points, and the number of the critical 
points can be zero, one, two, and three. 

Before presenting our numerical results, we shall give an intuitive account of the mechanism which causes the phase 
diagram with multiple critical points. We depict a schematic sketch in Fig. [3] which is useful to explain what would 
be anticipated in advance. There is a terminal point of the first-order phase boundary, that is the QCD critical point 
located at D, as long as the coupling ratio H/G is small and hence the neutrality constraint is not significant. When 
the coupling ratio H/G becomes substantially large, the critical line of the first-order phase transition is breached by 
the induced /i e in the CSC phase, and there appear two more edges in the critical line of the first-order transition; 
the new critical points are denoted as E and F. Here we notice that there exist three critical points as a whole. So to 
speak, the first-order line attached to F is a "survivor" transition that surpasses the diquark effect at sufficiently low 
temperature where the first-order chiral transition remains strong. The first-order line D-E is, on the other hand, to 
be regarded as a "remnant" where the diquark condensate almost melts and would hardly affect the chiral transition. 
For a larger H/G, the survivor may be gone and then the two critical points D and E arc left on the phase diagram. 
If we further increase H/G, all the first-order transitions of chiral restoration and all the critical points are washed 
away eventually. 

In the subsequent subsections we shall present numerical results and see what is described above is actually the 
case. For convenience we shall adopt the same notations as those in Ref. [|| to distinguish the different regions on 
the T-fi phase diagram; NG, CSC, COE, and NOR refer to the hadronic (Nambu-Goldstone) phase with a ^ and 
A = 0, the color-superconducting phase with A ^ and a — 0, the coexisting phase with a ^ and A / 0, and the 
normal phase with a = A = 0, respectively, though they have exact meaning only in the chiral limit. In fact, as seen 
from our results such as Figs. [5] (b), [7] (b), H and so on, M stays 10 ~ 100 MeV even in CSC but near COE. 

A. The Case of Intermediate Diquark Coupling 

We first consider the standard ratio H/G — 0.75, which is usually referred to as the "intermediate" diquark coupling 
strength. Two phase diagrams with and without the charge-neutrality constraint are presented in Fig. F|J henceforth, 
the thick solid curve, the dashed curve, and the thin solid curve in the T-/i plane stand for the critical lines of first- 
order phase transition, smooth crossover, and second-order phase transition, respectively. Since the magnitude of /xs 
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FIG. 4: The phase diagrams for the intermediate coupling constant H/G = 0.75 including the diquark condensate. The left (a) 
(right (b)) panel corresponds to the case without (with) enforcing the charge neutrality. In the left figure the unstable region 
is indicated by the dotted curve (see the discussion in Sec. IIII D|) . 
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FIG. 5: M, A, and \i e as functions of the chemical potential in (a) and the temperature in (b) for H/G = 0.75 under the 
condition of electric-charge neutrality. 



is very small [2l| , as we have mentioned, we ignore it in the numerical calculation. We have checked that nonzero /Js 
has only a minor effect on our main results. 

Both panels in Fig. |4] show that all the four phases, i.e., NG, CSC, COE, and NOR are realized in the phase 
diagram. The upper triple point at which the NG, CSC, and NOR phases encounter happens to be the critical point 
D for the chiral-to-CSC phase transition. 

Figure [4] (a) shows that the chiral phase transition keeps of first order at low temperatures even with the emergence 
of a COE phase when the charge-neutrality condition is not imposed. We remark that similar phase diagrams were 
obtained in the previous works [1, [U, [2f| [26| without imposing the charge-neutrality constraint in the NJL model, 
the instanton-based models, and the random matrix model. 
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In contrast to Fig. [4] (a), one finds an unconventional phase structure in Fig. [4] (b) where the charge- neutrality 
constraint is imposed: The would-be critical line for the first-order transition is terminated at two points E and F, 
between which the phase transition becomes crossover from COE to CSC. Hence there exist three distinct critical 
points, D, E, and F, as a whole. Although a possible phase structure with two critical points was suggested previously @, 
|8J, the present work is the first model study that shows possible existence of three critical points in the QCD phase 
diagram. 

In addition to the appearance of the three-critical-point structure, the following points are notable in Fig.F4](b): 

i) The critical line for the first-order chiral phase transition in the low-temperature region is shifted considerably 
toward higher quark chemical potential. Furthermore, the critical point D is somewhat shifted to lower T and higher 
fi; (T,n) = (42 MeV, 330 MeV) -> (33 MeV, 341 MeV). These are to be understood as the same effect of fi e as we 
mentioned in the explanation of Fig. [2] 

ii) The coexisting region is widened both in the T and fj, directions. The lower triple point of NG, COE, and 
CSC, labeled as E' in Fig. H (a) and E in (b), moves toward higher T and /x; (T, /i) = (5 MeV, 333 MeV) -> 
(25 MeV, 343 MeV). 

It is not difficult to understand the emergence of a smooth crossover for the chiral phase transition in the low- 
temperature region of the phase diagram when the diquark pairing is taken into account. The positive [i e first 
considerably reduces the discontinuity in the constituent quark mass at the first-order transition, and then the 
competition between the chiral and diquark condensates leads to the complete disappearance of the discontinuity. 
This situation is similar to what happens with the vector interaction as described in Ref. @; in Ref. Q, though 
the charge-neutrality condition is not imposed, a similar phase diagram to ours is obtained and the enlargement of 
the coexisting phase is attributed to the enhanced competition between the chiral and diquark condensations by the 
vector interaction. 

Now let us discuss the mechanism for realizing the three-critical-point structure shown in Fig. 0J For this sake, we 
also show, in Fig. [5j the /i and T dependence of M, A, and /i e in the upper and lower panels, respectively. We see 
the nature of the chiral phase transition from low to high temperatures in order. 

1) In the low-temperature region below the critical point F, we have a first-order transition from the COE to CSC 
phase, as shown in Fig. [4] (b). This feature is also clearly exhibited in the first panel of Fig. [5] (a) and the fourth 
panel in Fig. [5] (b) both of which show discontinuous jumps in the physical quantities. This result can be interpreted 
as follows; in this low-temperature region, the chiral condensate even in the COE phase has a rather large value and 
dominates over the diquark condensate at present diquark coupling H/G = 0.75. We should notice that the diquark 
condensate smears the Fermi surface like at finite temperature, which in turn tends to make the chiral transition 
weak [8]. In short, the effect of the diquark condensate is not yet strong enough in this temperature region to convert 
the first-order chiral phase transition to a smooth crossover. We identify this remaining first-order transition a survivor 
of the chiral transition which should be first order without the diquark condensate. 

2) Figure 0Jb) shows that at intermediate temperatures between E and F, the chiral transition becomes smooth. The 
underlying mechanism for this may be understood from an unusual temperature dependence of the diquark condensate 
in the relevant chemical potential region. In fact, as shown in the third panel of Fig. [5] (b), the diquark condensate 
increases as T is raised for a fixed fi. This is due to the combined effects of the charge-neutrality constraint and the 
temperature [2lj : The neutrality constraint causes mismatched Fermi spheres for the quarks involved in the pairing, 
which unfavors the pairing, especially at small temperatures. At larger temperature, the mismatched Fermi surfaces 
are smeared enough to allow for the significant number of the quarks involved in the pairing, and hence the diquark 
condensate can develop. Besides, the chiral condensate decreases in a monotonic way with increasing temperature, 
and a smaller quark mass favors the formation of the diquark condensate because of a larger Fermi surface at a fixed 
fi. Thus the effect of the diquark condensate may overwhelm the chiral condensate at intermediate temperatures, 
which suppresses discontinuity in the chiral phase transition which turns to crossover. 

3) In the still higher temperature region, we find that the phase transition comes back to first order from NG to 
CSC, which is between D and E in Fig. 2] (b). These features can be understood as follows. First, we note that D 
is attached to the critical line where A melts, and the diquark condensate decreases along D-E with the increasing 
temperature, as shown in the second panel of Fig. F|] (b). It means that the chiral condensate dominates over the 
diquark pairing in this region. Thus the original feature of the chiral transition without the diquark condensation is 
intact there, and COE does not appear. We see that the chiral restoration remains of first order as shown in the third 
panel of Fig. [5] (a) and the second panel of Fig. \5\ (b). In short, the first-order transition in this region is a remnant 
of the chiral restoration existing without the effect of the diquark pairing. 

4) For even higher temperatures than D, the chiral transition is a crossover and the condensates only show smooth 
behavior as seen in the fourth panel of Fig. [5] (a) and the first panel in Fig. [5] (b) where the second-order nature of 
the CSC transition is also exhibited. 
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FIG. 6: The phase diagrams for the strong diquark couplings H/G = 0.8 (a) and H/G = 0.875 (b) under electric-charge 
neutrality. The unstable region is indicated by the dotted curve as previously. 
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FIG. 7: M, A, and fj, e as functions of the chemical potential in (a) and the temperature in (b) for H/G = 0.8 under the 
condition of electric-charge neutrality. 



B. The Case of Strong Diquark Coupling 

We next choose H/G = 0.8 and H/G = 0.875 as examples for the "strong" coupling case, and present the corre- 
sponding phase diagrams in Fig.[6j One can see that, with increasing diquark coupling constant, the three-critical-point 
structure is first replaced by the two-critical-point structure and eventually the whole chiral phase transition becomes 
crossover. The phase boundary of the CSC phase is always a critical line of second-order transitions, which would be 
altered by gauge field fluctuations [27l ]. 

In contrast to Fig. [4] (b), Fig. [6] shows that, as H/G increases, the critical lines for the chiral and diquark phase 
transitions shift toward the lower chemical potential direction and the COE region is relatively enlarged. This is a 
natural result since the enhanced diquark condensate more strongly suppresses the chiral condensate. It is interesting 
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FIG. 9: The phase diagrams for the weak diquark couplings H/G = 0.7 and H/G = 0.6875 under electric-charge neutrality. 
The boundary of the unstable region is indicated by the dotted curve as previously. 




that the enhanced diquark condensate first converts the "survivor" line existing at intermediate coupling into crossover 
as noticed from Fig. [5] (a) and then melts the "remnant" line down to crossover as shown in Fig. [5] (b). This means 
that in the COE region the diquark condensate plays a dominant role over the chiral condensate and the original 
nature of the chiral transition fails to survive at low temperature. We remark that the two-critical-point structure in 
Fig. [5] (a) is topologically the same as what was reported in Ref. [1| , which may further illustrate an analogy between 
the positive fj, e and the repulsive vector interaction. 

The dependence of M, A, and \i e on \x (or T) for different fixed T (or fi) with H/G = 0.8 is presented in Fig. [7] (a) 
(or (b)). All the figures show that all the physical quantities behave in accord with the above picture. 

Figure [8] shows M, A, and /x e as functions of T at four fixed chemical potentials at H/G — 0.875. We find that both 
the remnant and the survivor are breached by the large diquark condensate and no first-order transition remains. 
We only notice that Fig. [5] shows again the unusual temperature dependence of the diquark condensate, especially 
in the COE region; the diquark condensate is an increasing function of T in the low-temperature region, and then it 
decreases in the higher temperature region. The mechanism for this peculiar behavior of the diquark condensate has 
been accounted for in the last subsection. It is worthwhile to note that the fourth panels in Figs. [7] (b) and [5] show 
that, if the diquark coupling is strong enough so that A ^ at T = 0, the nonmonotonous behavior of the gap energy 
exists no more, which only decreases with increasing temperature. 



C. The Case of Weak Diquark Coupling 

In this subsection, we present and discuss the phase diagrams for "weak" diquark coupling choosing H/G = 0.7 and 
0.6875. As we will see, weakening the diquark coupling can lead to unexpected complication in the phase diagram owing 
to the interplay between the chiral and diquark correlations, which is enhanced by the charge-neutrality constraint. 

Figure [9] (a) and (b) are the phase diagrams for H/G — 0.7 and 0.6875, respectively. Figure [9] (a) shows that the 
smooth crossover line E-F seen in Fig.|4](b) shrinks in a way that the critical point F meets E at E' on Fig. [9] (a). Then, 
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the thrcc-critical-point structure is replaced by the usual one-critical-point structure. This is because the magnitude 
of the diquark condensate is too small to turn the chiral transition into crossover. 

Let us briefly discuss a notable point that, for H/G = 0.6875, the critical point D and the triple point denoted by 
D' can become separate. It seems that, as long as the diquark coupling is larger than a certain critical value, there is 
a mechanism to make D and D' coincide. A similar observation is reported in the 2 + 1 flavor case too [28[. If H/G 
is lowered further, our numerical calculation results in the shrinking COE and the enhancing (low-T) NOR regions. 

Although we can recognize an intriguing structure in Fig. [9] (b) in comparison to (a) - the emergence of a NOR 
"island" surrounded by the CSC phase at low T - we will not take a close look into this region. This is because the 
newly arising NOR phase and the triple point G in Fig. [9] (b) are located deeply in the unstable region, the boundary 
of which is indicated by the dotted curve. We shall address the instability problem in the next subsection. 



D. Chromomagnetic Instability 

It has been known that the homogeneous CSC state could suffer from instability [2^, [3(| H3, HU, HH when the Fermi 
surface mismatch grows comparable to the pairing gap. The instability occurs in various channels simultaneously [34| . 
among which the transverse gluon field triggers the chromomagnetic instability. We can theoretically perceive the 
instability by the negative Meissner screening mass squared (or pure imaginary Meissner mass). In Figs. [4] (b), El an d 
IH1 we have located the unstable region given by the condition that either of eight gluons has negative Meissner mass 
squared. In what follows let us first explain how we calculate the Meissner mass. 

In the two-flavor case the analytical expressions for the Meissner mass are available in the hard dense loop (HDL) 
approximation [2^ . |35| . In our evaluation we have utilized them with our model parameters substituted to the 
formulae. We draw the dotted curves in Figs. HI (b), [H and [S] in this way. At the same time we have carried out 
the brute-force computation by means of the potential curvature with respect to the gauge field source [||| . This 
calculation is necessary for the numerical check of the consistency between the HDL in Ref. [35j and the approximation 
made in this work. We have in fact confirmed the consistency in the case that we keep using a sharp three-momentum 
cutoff A = 651 MeV. It should be mentioned that how to renormalize the unphysical cutoff dependence is not yet 
known in a field-theory manner [30]. Our prescription is motivated in accord with the NJL model treatment, while a 
larger A would make the unstable regions smaller. In this sense, even though the dotted curves in Figs. HI (b) , El [9] 
might move downward, the cutoff dependence would not affect the stable regions above them. 

We show the Meissner masses squared for the fourth and eighth gluons at H/G = 0.75 as functions of \x in Fig. [TO] 
We note that the Meissner mass squared changes rapidly from NG to COE leading to the chromomagnetic instability 
in the high-/z and low-T region of COE and CSC as indicated by the dotted curve in Fig. [4] (b). 

Figures [3] (b) and El (a) tell us that two critical points associated with the D-E line are outside the unstable region. 
However, the instability analysis solely can not tell the fate of another critical point F in reality. In fact, there are two 
possibilities to accommodate three critical points; one is that the unstable region may shrink to lower temperature if 
the cutoff is taken large and the other is that some other H/G value may push F up above the unstable region. The 
latter case can been confirmed in the current model study when H/G is located in a narrow region near but below 
0.75. To settle the robustness of F, we have to go into the identification of the true ground state inside the unstable 
region. 

There are some attempts to overcome the chromomagnetic instabi lity . The instability has a favorite direction leading 
to the (colored) Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) state [34|, Ha, [37| (or spontaneous current generation [38l . [39l |. 
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which shares the same mathematical structure with the colored plane- wave LOFF description [40|). We note that the 
gluonic phase [4l[ with one gluon condensate can translate into the plane- wave LOFF state, while a more stable gluonic 
phase with multiple gluon condensates has a chance to surpass the LOFF state 42]. The stable-unstable boundary 
we drew on the phase diagram can be regarded as a second-order phase transition line from the homogeneous CSC 
phase to one of these inhomogeneous states. Therefore, even though F in Fig. [4] (b) is overridden by the instability or 
the second-order phase transition, the physical consequence is very similar to the three-critical point structure, that 
is, one of the critical points is replaced by the critical line. Interestingly enough, then, Fig. [S] (a) is indistinguishable 
from Fig. F|J meaning that the three-critical point (or line) interpretation is in effect elongated for a wider range of 
H/G due to the chromomagnetic instability. 

E. Discussions 

In the previous subsections, we have given a detailed account for four possible scenarios for the critical point 
structure with the influence of the diquark condensate under the electric-charge neutrality. For the model parameters 
adopted here, numerical calculations suggest that the three-critical-point structure of the phase diagram appears in 
the range 0.735 < H/G < 0.767, while the two-critical-point structure in the range 0.767 < H/G < 0.82. For the 
stronger coupling region H/G > 0.82, the chiral phase transition remains crossover, while for the weaker coupling 
region H/G < 0.735, only one-critical-point structure appears. 

We have checked that the region with multi-critical-point structure is slightly modified quantitatively when nonzero 
fis is taken into account; the range 0.745 < H/G < 0.772 corresponds to the three-critical-point structure and 
0.772 < H/G < 0.80 to the two-critical-point one. Thus we can conclude that our main results on the novel phase 
structures are not altered even when the color-charge neutrality is fully incorporated although the parameter range 
leading to the multi-critical-point structure may be somewhat narrowed. 

It should be noted that the multi-critical-point structure was not pointed out in the previous works such as Ref. [28j 
where the electric neutrality and the dynamical quark mass were simultaneously considered for 2 + 1 flavor quark 
matter. Since the strange quark mass is still as large as the quark chemical potential near the chiral phase transition, 
the u and d quarks should play the dominant role in the vicinity of this region. Then, what causes the difference in 
the phase diagram? We have found that the appearance of the multi-critical-point structure is sensitive to the value of 
the constituent quark mass in the vacuum. In the present study the vacuum constituent quark mass for light quarks 
is 325.5 MeV, while it is 400 MeV in Ref. f28j . Notice that the larger the constituent quark mass is, the smaller the 
Fermi sphere becomes, disfavoring the diquark pairing for a given fi. Therefore, in Ref. [28| . the strong first-order 
chiral transition at low temperature is hardly affected by the diquark even with fj, e , which could make the chiral 
transition smooth. However, in turn, it suggests the possibility that the multi-critical-point structure may appear 
even with the large constituent quark mass adopted in Ref. [28j | if the repulsive vector interaction is also included, 
which is rather close to the realistic situation. 



IV. CONCLUSIONS AND OUTLOOK 

In this paper, we have explored the effect of electric-charge neutrality in (3 equilibrium on the chiral phase transition 
with and without considering the diquark condensate within a simple two-flavor NJL model. 

We first disclosed the similarity of the roles on the chiral phase transition played by the positive electric chemical 
potential and the repulsive vector interaction as follows: 1) Both of them make the chiral critical line at low temper- 
ature shift toward larger quark chemical potential and effectively weaken the first-order chiral phase transition. 2) 
Both of them effectively enhance the competition between the chiral and diquark condensates. 

In some model parameter region including the diquark condensate, the combination of these two properties can 
result in the emergence of the two-critical-point structure for the chiral phase transition in the electric-charge neutral 
case. It is noticeable that the magnitude of the repulsive vector interaction in Ref. [8] was an assumption although 
the existence itself of such an interaction has a generic foundation; on the other hand, the positive electric chemical 
potential in this paper is self-consistently determined by the physical requirement that the bulk matter must be under 
the charge-neutrality condition. 

Besides two nontrivial effects mentioned above, our investigation also showed a quite unconventional three-critical- 
point structure for the chiral phase transition. This result is directly associated with the abnormal behavior that 
the diquark energy gap can increase with increasing temperature in the coexisting region for a certain range of T-fi 
and H/G, which is possible under the tangled influences of electric-charge neutrality and the interplay between the 
chiral and diquark condensates. To our knowledge, so far, this is the first case which concretely demonstrates the 
three-critical-point structure for the chiral phase transition. Although this result with three critical points could be 
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sensitive to the model parameter choice, it definitely gives a hint that there may exist a possibility that the complexity 
in QCD allows for more complicated phase boundaries associated with the chiral phase transition. Moreover we have 
studied the chromomagnetic instability and found that the three-critical-point scenario may be taken over by the 
structure with two critical points and one critical line. Our result also has a meaningful implication for the study of 
phase transitions leading to a coexisting phase involving superconductivity in condensed matter physics when external 
conditions are enforced to the system like the neutrality constraint causing the mismatched Fermi surfaces. 

In view of the common effects on the chiral phase transition by the positive y, e and the repulsive vector interaction 
Gy, it is natural to expect that the influences mentioned above will be greatly enhanced when these two aspects 
are considered simultaneously. It is expected that the three-critical-point structure may extend to a weaker coupling 
region in the presence of the vector interaction, although the number of the critical points might change in a different 
way from that depicted in Fig. [3] when the vector coupling is increased. 

In addition, the confinement-deconfinement phase transition at finite temperature and density is still poorly under- 
stood inside the CSC phase, that is because of the lack of confinement in the NJL-type model. Recently the PNJL 
was proposed 1-43 . I3| and had been extensively used to study the thermal properties and the phase transition of 

qcd us, m snuiss ee h, m, m mi. It is interesting to investigate whether the results obtained in this paper 

are also true when including the Polyakov loop dynamics. Practically, the PNJL model tends to push the location of 
the critical point up toward the higher temperature region. That means, a detailed structure observed in this work 
could be magnified by the Polyakov loop. More details on these two issues will be reported in our future publica- 
tion [l5j . and in fact, we have already verified that the multi-critical-point structure still exists when the Polyakov 
loop effect is taken into account. 
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